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Abstract. We investigate the bias dependence of the tunneling conductance 
between a spin-polarized (SP) scanning tunneling microscope (STM) tip and the 
surface conduction states of a normal metal with a Kondo adatom. Quantum 
interference between tip-host metal and tip-adatom-host metal conduction paths 
is studied in the full range of the Fano parameter q. The spin-polarized STM 
gives rise to a splitting of the Kondo peak and asymmetry in the zero-bias 
anomaly depending on the lateral tip-adatom distance. For increasing lateral 
distances, the Kondo peak-splitting shows a strong suppression and the spin- 
polarized conductance exhibits the standard Fano-Kondo profile. 
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1. Introduction 

The Kondo effect is an antiferromagnetic screening of a localized magnetic moment 
by the host metallic electrons below a characteristic Kondo temperature Tk, which 
results in the appearance of an additional peak in the system's density of states pinned 
to the Fermi energy. Being first observed during studies of transport properties of 
bulk diluted magnetic alloys pQ, the Kondo effect was later on shown to affect the 
conductance of single quantum dots (QDs) [H13], arrays of QDs [U [5] and possibly 
quantum point contacts [7] . In the emerging field of spintronics [8] , the coupling 
of a single QD to ferromagnetic leads can shed more light on fundamental aspects of 
the Kondo physics and provide a basis for a design of novel spintronic composants 

la uni nn m m m ee ces on . 

The Scanning Tunneling Microscope (STM) is widely used for investigation of the 
Kondo effect at surfaces of normal metals with adsorbed magnetic impurities (Kondo 
adatoms) [18l[l9l[2Ql[ni[22[23[2ll[25l[26]. The possible examples are individual Co 
atoms at Au(lll), Cu(100) or Cu(lll) surfaces QH Q3B |13 HD , and cobalt carbonyl 
Co(CO)„ complexes or manganese phthalocyanine (MnPc) molecules on top of Pb 
islands [221123] , 

In the present work we focus on the Kondo regime for a system containing a 
ferromagnetic STM tip and a single Kondo adatom on a metallic surface (figure [T]) . 
In this system the conductance becomes spin dependent, due to the ferromagnetic 
(FM) tip, which leads to a Zeeman splitting in the Kondo adatom density of states 
(DOS). Particular attention is payed to the interplay between Kondo effect, quantum 
interference of two possible tunneling paths (tip-adatom-host and tip-host) and the 
ferromagnetism of the tip. 

Our study differs from previous work in the same field [27j mainly in the three 
aspects described below. 

First. We present new results for the tunneling conductance as a function 
of the bias for different lateral tip-adatom separations R in the Kondo regime. 
This is an actual task indeed, as in experiments with unpolarized systems, R was 
shown to strongly affect the conductance [HI [19] . For ferromagnetic tips, however, 
correspondent results are lacking. 

Second. We discuss the small, large and intermediate cases for the Fano parameter 
q [28l [29l [30l [3T], [32] . We show that the three cases above give different conductance 
patterns. Experimentally the large and small q limits are relevant and have been 
realized for unpolarized conduction bands. For example, in reference [21] it was 
shown that q alternates from intermediate (q « 1.13) to small (q w 0.18) values 
by changing the Cu crystal surface (Cu(100) and Cu(lll)), in which a Co adatom 
is deposited. The conductance profile for small q limit can also be experimentally 
observed in quantum corrals [19] . Alternatively, manipulating the molecular structure 
of the magnetic impurity on the Cu(100) surface, it is possible to switch between the 
intermediate and large q limits (22[ [23] . 

Third. Instead of using a Lorentzian approximation for the Kondo peak, we 
apply the well established Doniach-Sunjic formula [331 134] . whose validity is supported 
by numerical calculations based on both numerical renormalization group [33|, 134] . 
quantum Monte Carlo simulations [35] and by the recent success in fitting experimental 
data [23] . 

Our results show that for a tip situated right above the Kondo adatom (R=0) 
an asymmetric zero-bias anomalies appear, which are revealed as resonances and anti- 
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Figure 1. Setup of the SP-STM and the host metal with the Kondo adatom. 
The interference between the channels t^n and t c exhibits a spin-polarized Fano- 
Kondo profile for the conductance. 

resonances in the conductance in the limits of large and small q, respectively. For 
q rs 1, the conductance demonstrates a pronounced plateau in the region of small biases 
(eV ~ 0). The increase of R leads to a suppression of the tip-induced adatom's Zeeman 
splitting, thus resulting in a conductance pattern that resembles experimental data 
for unpolarized systems |18j . Finally, we verify that the Kondo peak splitting strongly 
depends on the asymmetry between the tip- adatom and adatom- host tunnelings. 

The paper is organized as follow. In section 1 we present the model adopted 
to describe the system under study. In section 2 we derive an expression for 
the spin-resolved conductance based on a perturbative expansion for the tunneling 
Hamiltonian. In section 3 we discuss the numerical results. Conclusions are present 
in section 4. 



2. The Model 

The description of the metallic host and its interaction with the adatom is performed 
within the framework of the single impurity Anderson model [36) with a half-filled 
noninteracting conduction band, 

Ha= ^2 ec\ a c Ea de + ^ e Aa d\d a + Un d ^n dl 

<7 a 

+ /^/"de( C U + ff. C .), (1) 

where all energies are measured from the Fermi level coinciding with the center of the 
band (ef = 0) and extend from — D to D. The operator 

c ea = Po 1/2 ^2cj:J(e ~ e k ), (2) 

k 

corresponds to a surface conduction state of the host metal with an energy-independent 
density of states per spin p . 

The first term in equation (Q} describes a two dimensional electron gas (2DEG) on 
the surface of a host metal. The second and the third correspond to the adatom, which 
is characterized by a single particle orbital energy Eda and Coulomb repulsion U. The 
last term, proportional to vT 7 , describes the coupling between the adatom and the host 
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metal conduction states, thus introducing the broadenings of the adatom's resonances 
at the energies and e^+U. In the Kondo regime (T <C T K , e da < e F , Sd a + U > £p, 
r <C \£da\, £da + U) an additional peak in the density of states having a half-width 
r_ff = ^bTk, (fcs is the Boltzmann constant and Tk the Kondo temperature) [37J 
appears exactly at the Fermi level. 

The total STM Hamiltonian reads 



TM 



Ha + H t i„ + H t 



where H t i P corresponds to free electrons in the tip, 



(3) 
(4) 



with the operators ap a describing the bulk conduction states, the bias is eV and 



Htun = J2 [*c4r-M#) + H C - 



(5) 



is the tunneling Hamiltonian that connects the tip with the host metal via the operator 



B a (R) 



The normalization factor in equation ^ reads 



N K = 



-1/2 



(6) 



(7) 



with tpz(R) ~ e being a wave function of the host conduction electron. 
The first term in equation ^ containing the operator 



(8) 



hybridizes the conduction states of the tip with the surface of the host at a displaced 
lateral position R from the adatom. 

The second term describes the tunneling between the tip and localized adatom's 
level characterized by a spin-dependent parameter (Fano factor) [351 [3S] 

q ° R = (nr Po y 1/2 (t° dR /t°) , (9) 

defined as a ratio between the couplings t a dR of the tip-adatom and t% of the tip-host 
metal. 

The Fano factor monitors the competition between the tunneling channels in the 
system. It vanishes with the increase of the tip-adatom lateral distance R, which can 
be modeled by an exponentially decaying function |39j 

qn = q R=0 e- k - R . (10) 

The decay of the Fano parameter with the increasing of the lateral tip-adatom 
distance was already experimentally explored in a system with a Co adatom on a Cu 
surface [21]. 

In the present work we consider a ferromagnetic tip with a spin-dependent density 
of states given by 

P? ip = Po [1 + <rPti P ] , (11) 
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where a = + or — for spins \ or j, respectively, P^p is a polarization degree of the 

T I 

tip. The inequality between spin-up and spin-down populations in the tip p tip > p tip 
introduces an asymmetry in the splitting of the zero-bias anomaly as we will see in 
Sec. IV. The density of states for the unpolarizcd tip po = ptipi^F) is assumed to be 
equal to the density of states of the host metal for simplicity. 



3. Methodology 



We calculate the tunneling conductance of the system treating the coupling between 
tip and host metal (H tun ) as a perturbation. Within a second order perturbation 
scheme, the formula for the conductance EE] reads 



G=(e 2 /h)J2 J T a (e,T,R,q R ) 



d_ 

de 



f(e - eV) 



ds, 



where 



(e,T,R,q R ) = T oa [l + \q R \ 2 } (f% p /po) (pIdos/po) , 



is an effective transmission coefficient with 

,2 



(12) 

(13) 
(14) 

(15) 



The local density of states (LDOS) appearing in equation (TT3|) is defined as 

l %((B a (R)\Bj(R))) E 

PLDOS ~ * i + Kf ' 

with ((B cr (R)\B'l(R))) £ being the retarded Green function thermally averaged over 
the eigenstates of the Hamiltonian ^ . 

Looking at equation (|15p. one sees that the LDOS depends on non-orthonormal 
fermionic operators 



where the spatial function F^(e) is given by 

f r( £ ) = £ n(R) 5 ( £ ~ £ k) = p Jo(k{s)R), 



(16) 



(17) 



with Jo being the 0-th order Bessel function. 

The operator ^ can be expressed in terms of fermionic operators orthonormal 
to c ea by introducing 



N 



with a normalization factor evaluated at the Fermi level, 



^Vo = < 1 - 



-1/2 



(18) 



(19) 
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This leads to the following expression for B„(R), 



B a (R) 



J dec ea 
NoNo 



Pi) 



where we introduce an operator 



r de, 



(20) 
(21) 



which describes a conduction state centered at the adatom site. 

As the Kondo effect occurs at low temperatures T <C Tk and the tip bias is 
usually much smaller then the bandwidth, eV <§; D, we evaluate equation (fT2)) at 
T = 0, thus resulting in 

G(eV, T -c Tk, R) — ^ G a max (pu p /po) (pIdos/po) 



where 



and 



G" 



(e 2 /h)T oa {l + \q R \ 2 } 



(22) 
(23) 



PIdos/po = [1 - Jo (k F R)] cos 2 % R + sin 2 6° H sin 2 % v 
+ 2 J (kpR) sin 5 a qR cos 5° R sin 5° v cos 5 a eV 
+ J 2 (Jbf-R) cos 2 S° R COS 2 tfJvr. (24) 

In the calculation of the above expression we used the Green function identities for 
zero temperature Anderson model [36] . The spin dependent Fano factor phase shift 
8°„ is defined as 

t & n5° R = \<f R \. (25) 

In equation (|24[) the terms proportional to sin 2 5° v and cos 2 6° v come from the 
direct tunneling paths tip-adatom-host and tip-host, respectively. The interference 
between then is given by the mixture term proportional to sin^y cos 6°y- 

The spin dependent phase shift 5° for the conduction states can be determined 
from the Doniach-Sunjic spectral density (33J, [34] 



pUeV) 



1 



3? 



iT 



K 



(eV + crA) 



iT K 



1 

ttT 



= sin 2 Sly, 



(26) 



where A gives the Kondo peak splitting, which is related to the adatom energy level 
as 



o A 
£da =e d -a—, 



(27) 



where e° d is the adatom level without magnetic field. According to "poor man's'' 
scaling [TT], such splitting can be estimated as [ID] 



A _ tip 



1 tip 



2tt 



^ \n(D/U) 

7T 



P tlp MD/U) 

Pti P exp(-2k F R) . 



(28) 
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Figure 2. Conductance G/Gmax as a function of the tip bias scaled in units of 
Kondo resonance half- width eV/Tjc, for qR—d = 10 at three different tip lateral 
positions with symmetric potential barriers. 

where T® ip = tt \ t^ R=0 \ 2 po gives the local coupling between the Kondo adatom and an 
unpolarized tip. We use the expression 

rSp = r°, lp [1 + aP tlp ] exp (-2k F R) , (29) 
to account for both spin and spatial dependencies of the tip-adatom coupling. 

4. Results 

For numerical analysis we adopt the following set of model parameters: e° d — — 0.9eV, 
r = 0.2eV, U = 2.9eV, D = 5.5eV, T K = 50K and k F = 0.189A" 1 [2SH7|. We 
consider the cases of large, small and intermediate Fano ratio values. For each of these 
cases we analyze the dependence of the conductance on tip-adatom lateral distance 
R. Additionally, for large q limit we consider how the conductance depends on the 
asymmetry between tip-adatom and adatom- host couplings. 

4--1- Large tip-adatom coupling (qn = o = 10) 

The large q limit has been achieved in experiments with STM tip by employing 
magnetic molecules as Kondo adatoms UH [23] ■ In this limit the host metal conduction 
electrons tunnel towards the tip preferably via the localized magnetic adatom state. 
For the tip situated right above the adatom (R = 0), the conductance reveals an 
asymmetric splitting of the zero bias anomaly (figure (HJ)), characterized by a pair of 
peaks at eV — —A and eV = A. Such asymmetry occurs due to the spin polarization 




Figure 3. Conductance G/G m ax at the origin as a function of the tip bias scaled 
in units of Kondo resonance half- width eV/Tx, for qji=o = 10 an d different 
potential barriers. 



of the tip (equation (jTTJ) ) , the higher peak corresponding to the majority spin up 
states, while the lower one to the minority spin down states. 

A similar asymmetric splitting of the Kondo peak was recently observed in a 
quantum dot system coupled to two ferromagnetic Ni electrodes [151 E] • The large 
q limit in the system we consider resembles the standard case of a single dot in 
between leads without lead-to- lead direct coupling (embedded geometry). A wealth 
of theoretical works predict the spin splitting of the Kondo peak in a quantum dot 
system coupled to two ferromagnetic leads [TU1 EJ [T5] • In those works this splitting 
was tuned via the relative angle between the left and the right lead magnetization 
[T5] . the leads polarization 10J and an external magnetic field [TT]. Here we show 
one alternative/additional way to tune the spin splitting, by changing the tip-adatom 
separation (laterally or vertically) [41] . 

Increasing the tip-adatom lateral distance, the Fano ratio decays according to 
equation (fTU)) and the Zeeman splitting of the Kondo peak quenches (see equation 
(|28p). These effects can be seen at figure ©, where the conductance is plotted for 
three different values of the tip-adatom lateral distance. About R = 6.65A, the two 
resonances merge into a single peak, thus resulting in the standard Kondo resonance 
profile. 

Not only the lateral tip-adatom separation can change the spin splitting, but also 
the vertical tip-adatom distance. To mimic this effect we change the ratio T^ ip /T, 
thus introducing an asymmetry between the tip-adatom (T^ ip ) and adatom-host (r) 
tunneling rates. When the tip becomes vertically closer to the adatom, we expect the 




increase of the coupling parameter Consequently, the splitting between the peaks 
becomes more apparent (see figure ([3])) which can be understood as consequence of the 
enhancement of the local tip magnetic field on the adatom due to the tip proximity. 



4-2. Small tip-adatom coupling (qn = o = 0.01 ^ 

In the small coupling limit the conductance curves (figure display dips instead 
of peaks observed in the large coupling limit (figures (J2j)-(j3j) ) . The appearance of the 
dips at eV = ±A is a consequence of a destructive quantum interference between 
the channels tip- host and tip-adatom- host. This can be easily seen from equations 
JM])-([23 for q -> and R=0 that the LDOS behaves as 

pIdos/po = cos 2 5: Vl (30) 

which is opposite to the case qR=o — > oo, where 

Pldos I Pa = sin2 Kv ( 31 ) 
As in the case of large coupling, the increase of the tip-adatom distance leads 
to a quenching of the anti-resonance splitting as it can be seen comparing the curves 
corresponding to different values of R at figure ((4]). At R — 6.65A, the asymmetric 
zero-bias anomaly for the dips disappears and only the standard single anti-resonance 
profile is recovered. A single dip structure in the conductance for small q is verified for 
a system composed of Co on Cu(lll) surface [TH1I2I]. It is valid to note that while the 
large q limit resembles the embedded geometry, the small q limit gives similar results 
to the T-shaped quantum dot (side-coupled geometry) [3l |42l [43] . 
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4-3. Intermediate tip-adatom coupling (qu=a = 1) 

The intermediate case for a normal tip is characterized by the well known Fano-Kondo 
line shape [181121] . However, the introduction of a ferromagnetic tip results in distinct 
Fano-Kondo profiles for each spin component. The spin up profile is shifted toward 
negative bias with an enhanced amplitude, while the spin down case moves in the 
opposite direction being reduced in amplitude. This is illustrated in the inset of figure 
iJSJl, where we show /G max and G^/G max - 

The superposition of the spin-dependent Fano-Kondo profiles gives rise to the 
appearance of a plateau around the Fermi level (eV = 0) in the total conductance 
(G^ + G^) for R = (see figure ([5])). For large R (R — 6.65A), this plateau vanishes 
due to the suppression of the adatom's Zeeman splitting and standard Fano-Kondo 
results for the conductance with nonmagnetic tips are recovered [T81 H?T] . 

5. Conclusions 

We derived a spin resolved tunneling conductance for a system of a spin polarized 
STM tip with a single Kondo adatom on the surface of a normal metallic host. The 
conductance dependence on the tip bias was investigated for different tip-adatom 
lateral distances in a wide range of the Fano parameter q, relevant for various 
experimental configurations. We demonstrated that the Fano parameter drastically 
affects the conductance pattern of the system. For large values of q, we observe an 
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asymmetric splitting of the Kondo resonance which is suppressed with an increase 
of the tip-adatom lateral distance. For small q, the behavior of the conductance is 
opposite to those observed in the large q regime- instead of a splitted Kondo peak, one 
observes an asymmetrically splitted Kondo dip. For the intermediate case we have 
shown that due to the splitting of the spin resolved conductances (G^ ^ G the 
total conductance exhibits a plateau in the region of small biases, which disappears 
for large enough values of a tip-adatom lateral distance. 
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